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Abstract 


A bivariate  failure  model  is  proposed  in  which  the  residual 
lifetime  of  one  component  is  dependent  on  the  working  status  of 
the  other.  General  properties  of  the  model  are  discussed,  and 
the  maximum  likelihood  estimates  of  the  parameters  are  found 
in  a bivariate  exponential-like  special  case. 


Keywords:  Bivariate  failure  model,  bivariate  exponential 

distribution,  maximum  likelihood  estimation. 


I failure  of  one  component  puts  more  (possibly  less)  strain  on  the 

I remaining  components,  for  example,  the  kidneys.  Section  two 

derives  properties  of  the  lifetimes,  including  their  joint  Laplace- 
I Stieltjes  transform. 

In  the  third  section  a bivariate  exponential- like  special 
[ case  is  considered.  In  this  example  maximum  likelihood  estimates 

; of  the  parameters  are  obtained  and  their  asymptotic  distribution 

[ studied.  This  model  is  compared  with  the  bivariate  exponential 

I models  of  Freund  [3]  and  Marshall  and  Olkin  [41. 

! 

I 

2.  MODEL  DEFINITION  AND  GENERAL  PROPERTIES 

Label  the  two  components  of  the  system  A and  B with  lifetimes 
S and  T respectively.  The  lifetimes  of  the  two  components  are  de- 
pendent, in  that  the  failure  of  one  component  affects  the  residual 
lifetime  of  the  other.  To  describe  S and  T,  let  X,  Y,  U,  V be 
non-negative  mutually  independent  random  variables  with  X and  Y 
absolutely  continuous.  Then  we  write 

S = min(X,Y)  + U • Ifx>  yi 

(2.1) 

T = min(X,Y)  + V • Yj 


We  will  obtain  an  expression  for  the  joint  survival  distribution 


2 


F(s,t)  = Pr{S>s,  T>t}.  Assuming  s<t  and  conditioning  on  the 
values  of  X and  Y,  we  obtain 


F(s,t) 


-OO  -CO 

QJ  0 


Pr{S  > s. 


T>tlx=x,  Y=y}dFj^(x)dFY(y)  . 


Now  partition  the  region  [0,«>)  x [0,°°)  into  12  subregions  based  on 
the  relative  sizes  of  x,y,s>t.  The  only  non~zero  contributions  are 
PrfV  ; t-xl  on  the  region  [s,t]  x [x,«>)  and  1 on  the  region  [t,«>)  x [t,«> 
Thus  the  following  can  be  established. 


Theorem  1 : 


F(s,t)  = 


F^{t)F^(t) 

> Fj^(s)Fy(s) 
Fj^(s)  Fy(s) 


F^(t-x)FY(x)dFj^{x) 


+ j Fy(s-y)Fj^(y)dFY(y) 


if  s < t 
if  s = t 
if  s > t 


The  joint 
be  f * (a ,b)  = 
determined  by 
represented  by 


Laplace-Stielt jes  transform  of  (S,T)  is  defined  to 
as-btp (ds,dt) , where  F(ds,dt)  is  the  measure 

O-*  0 

the  survival  function  F(s,t).  This  measure  can  be 


F (ds ,dt) 


FY(s)dF^(t-s)dF^(s)  if  s < t 

• FY(s)p^dFj^(s)  + Fj^(s)PydF„(s)  if  s = t 

Fj^(t)dFu(s-t)dFY(t)  if  s>t, 


(2.2) 


where  = Pr{V=0)  and  Py  = Pr{U=0}.  This  expression  is  used  to 
evaluate  f*{a,b)  and  we  get 
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Theorem  2 : 


f*(a,b)  = f*(b) 
V 


e (s)dF  (s)  + f*  (a) 

) Y X U 


e ®F  (s)dF  (s)  . 

0 X Y 


These  integrals  cannot  be  evaluated  in  general,  but  can  be  for 
certain  important  special  cases.  In  particular  we  have 


Corollary  3:  If  Y has  an  exponential  distribution  with  parameter 
3,  then 

f*(a,b)  = f* (b) f* (a+b+B)  + (f * (a) B/a+b+6)  [1-f * (a+b+B) ] . 
Moments  can  be  calculated  by  differention  using  this  expression. 


3.  ESTIMATION  IN  A SPECIAL  CASE 
As  a special  case  suppose  that  X and  Y are  exponentially 

“ Cl  * t 

distributed  with  parameters  at  and  B respectively,  Pr{U  ’ t}  = qe 
S ^ t 

Pr{  V > t } = qe  , where  a,B,  ot',  6*'0,  O^q^l,  t_0.  Freund  ' s 
[3]  model  corresponds  to  the  case  q=  1.  The  parameter  q allows 
for  simultaneous  failure  of  the  components,  since  Pr{S=T}  = 1-q  = p. 
Using  results  from  section  2,  properties  (including  moments)  can 
be  derived.  In  particular. 


F(s,t) 


^-(a  + 6)t  ^ qae 


-e't  r 


ni+B-B 


-(ot+8-B  )s  -(a+B-6  )t 
e -e 


G 


- (u+B) s 


-(oi+B)s  qBe  ™ ® F - (nt  + B-oi ')  t - (n  + B-n  ' ) s H 

a+H-it  ’ L?  J 


a+B-u 


if  s < t 
if  s = t 
i f s t 


(3.1) 


Marshall  and  Olkin  [4]  defined  a bivariate  lack  of  memory 
property  by  PrlS  > s+A,  T > t+A}  = Pr{S>s,  T>s}'Pr{S>A,  T>A} 
for  all  s,  t,  A 0.  The  survival  distribution  (3.1)  possesses 
this  property  and  has  mixtures  of  exponential  distributions  as 
marginals . 

The  measure  determined  by  (3.1)  however  is  not  absolutely 

continuous.  If  we  let  represent  Lebesgue  measure  on  , then 

the  measure  defined  by  u (A)  = U2  (A)  + ij^{x:(x,x)  c A}  is  a suitable 

dominating  measure  for  maximum  likelihood  estimation.  This  is 

the  same  dominating  measure  used  by  Bhattacharyya  and  Johnson  [1]. 

Consider  a sample  of  N independent  observations  on  (S,T), 

{ (S  , t, ),...,  (S„, t„)  } . The  following  notation  simplifies  the 
1 1 N N 

likelihood  function . Let  A,  = {{S.,t.)|s.<t.},  A_  = {(s.,t.)|s.>t.} 

1 i i'i  i 2 i i 'i  1 

and  A^  = {(s^,t^)|s^  = t^}.  further  let  N^=  #A^,  N2  = #A2,  N2  = #A2, 
S=ys.,  S,  = ^s.,  R=^s.,  T =Jt.,  T^  = ^t..  Here  #A  is  the  number 

^a1  a1  1a1  ^a^ 


of  items  in  A.  Then  L,  the  likelihood  function,  can  be  expressed 


L = (aS')  Nl(i.p)Ni+N2pN3 

( a+B)  ( Bu  ' ) ^2exp  [ - (n+B)  (S j^+R+T2 ) “ 


a'  (S2-T2)-6'  (Tj^-S^)  ] . 


(3.2) 


Solving  for  the  maximum  likelihood  estimates,  one  obtains 


Theorem  4; 

(i)  If  Nj^=N2=0,  then  p=l  and  a,  6,  a',  6'  cannot  be  estimated, 


If  N^=0,  but  then  p=N2/N,a  =0,6  =N/ (S^-i-R+T2 ) , 

a'=N2/ (S2-T2),  and  6' cannot  be  estimated. 

If  Nj^^O,  but  N2  = 0 then  p=N2/N,  6 = 0,  a=N/(S^+R+T2)  / 

0 "=N^/ (T^-S^)  , and  u' cannot  be  estimated. 

If  N^?^0  and  N?«0  then 

a=(N/S^+R+T2) (N^/N^+N2) 

6=  (N/Sj^+R+T2)  (N2/N^+N2) 

a^=N2/(S2-T2) 

6'=Nj^/(T^-S^) 

P=N3/N. 

It  is  easy  to  obtain  the  biases  in  these  estimators.  Using 
a Lehmann,  Scheff4  partitioning  operation  (cf.  Zacks  [5],  p.  50) 
it  can  be  shown  that 

Theorem  5:  The  vector  (N^,  N2,  Sj^+R+T2,  S2-T2,  T^-S^)  , (and  thus 

the  vector  of  maximum  likelihood  estimators)  is  a minimal  sufficient 

statistic  of  the  sample  { (s. , t, ),..., (s„ , t„) } . 

i 1 N N 

To  obtain  the  asymptotic  distribution  of  the  maximum  like- 
lihood estimates  we  need  to  restrict  the  parameter  space  as  follows. 
Let  fi*  = {(a,6,a',6',p)|0<tj^<a<M^,...,0<t3<p<M3<l}.  On  this 
space  the  regularity  conditions  presented  by  Chanda  [2,p.  56]  are 
satisfied.  Thus  we  have 

Theorem  6: 

(i)  (a,  6,  a',  B'tP) — ^(a,  8,  a',  p)  as 


(ii) 

(iii) 

(iv) 


with  probability 


I 


< 

i 


I 

! 


i5 


i 

6 

• /V  - A.  /N 

(ii)  N'(a-a,  j'-ti",  B'-B',  p-p)  is  asymptotically 

distributed  as  multivariate  normal  with  mean  0 and  covariance 
matrix  'l , where 


-1 


a+B-p6 
=;(a+B)  ' 


(a+B)  ^ 


g+B-pu 


6 (a+B) 


(1-p) 


(a  ')  (a+6) 


(1-p) g 


(B')  ' (a+B)^ 


li 


P (1-p) 


The  parameters  can  also  be  estimated  by  maximum  likelihood  for  a j 

I 

censored  sample. 

V7e  see  that  this  model  generalizes  Freund's  model  to  include 
simultaneous  failure  of  both  components.  It  differs  from  the  model 
of  Marshall  and  Olkin  in  that  the  residual  lifetime  of  one  component 
is  not  independent  of  the  status  of  the  other  component.  We  feel 
that  these  features  will  aid  in  the  application  of  the  model. 
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